Coexistence of spanning clusters in directed percolation 
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The probability distribution for the number of top to bottom spanning clusters in Directed per- 
colation in two and three dimensions appears to be universal and is of the form P(n) ~ exp(— an 2 ). 
We argue that a is a new critical quantity vanishing at the upper critical dimension. The probability 
distribution of the individual masses of the spanning clusters is found to have a Pearson distribution 
with a lower cutoff. Various properties of the clusters are reported. 
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In a large variety of equilibrium or noncquilibrium 
critical systems spatial or space-time anisotropy plays a 
crucial role. Many such athermal cases exhibit a con- 
tinuous transition in the geometric properties like the 
sizes or shapes of clusters, and the critical behavior is 
best described as belonging to the universality class of 
directed percolation (DP) [Q. Examples abound Q as 
fluid flow through porous medium in an external field, 
threshold dynamics in random systems, forest-fire or epi- 
demic models, reaction diffusion systems, self organized 
criticality, resistance(or insulator)-diode network, dam- 
age spreading, Reggeon field theory etc. Naturally the 
detailed properties right at the percolation threshold p c 
arc influenced by the nature of the spanning clusters. 
But, then, how many spanning clusters does one expect 
at the threshold? 

Let us consider a DP problem on a square or a cu- 
bic lattice where each lattice site is occupied indepen- 
dently with a probability p. We want percolation that 
involves a top to bottom spanning (TB-spanning) clus- 
ter. Above the percolation threshold, there is only one 
infinite cluster occupying a finite density of sites of the 
lattice. However, at the threshold, a spanning cluster is a 
zero-density self-affine fractal (requiring direction depen- 
dent length rescaling). The lattice can then in principle 
support many such TB-spanning clusters. 

Although renormalization group arguments have led 
to a coherent understanding of long distance (and time) 
properties of DP, especially the universal exponents, de- 
tailed questions like the number of spanning clusters at 
the critical point still remains unanswered. The known 
exponents do tell us that there is a problem in high 
enough dimensions due to violations of the hyperscal- 
ing relation and this can be resolved |3],[iJ if there is an 
infinite number of spanning clusters at p — p c . This ar- 
gument is not applicable for dimensions below the upper 
critical dimension (d c = 5) and therefore the question 
stays whether there will be a finite number of spanning 
clusters or a finite density or just one. 

In a number of recent studies p|-|lT| , it has been found 
that there exists more than one TB-spanning cluster at 



the percolation threshold for d < 6 for ordinary (i.e., 
undirected or isotropic) percolation (OP) Jl2]. For two 
dimensions, it has been shown analytically that there ex- 
ists more than one TB-spanning cluster at the percola- 
tion threshold and the probability P(n) of the existence 
of n such clusters follows a behavior like exp(— an 1 ) with 
7 = 2 for large n [pi. There is some controversy about 
the behavior of P(n) in higher dimensions, though recent 
results indicate that 7 is independent of dimension and 
remains equal to 2 while a decreases with dimension [9| . 
The individual masses of the spanning clusters have been 
shown to follow the scaling behavior |l3| of the unique 
spanning cluster in all dimensions 2 to 5. It has been con- 
jectured that in two dimensions, the ratio of the masses 
of the largest and the second largest spanning cluster is 
in all probability universal |Tcj . 

In sharp contrast to ordinary percolation p4[| , very few 
exact or rigorous results are available (to the best of our 
knowledge) for DP (except mean field results). In order 
to develop an analytical understanding, attempts have re- 
cently been made to bridge the gap by studying crossover 
from OP to DP . Hence the necessity to know better 
the difference or similarity between OP and DP in fea- 
tures that are universal but yet to be resolved in the RG 
framework. 

Here we have studied coexisting TB-spanning clusters 
by simulating DP on square and cubic lattices with heli- 
cal boundary conditions in the transverse directions. The 
clusters are identified by using the Hoshen Kopelman al- 
gorithm, respecting the directedness of the paths. For 
any finite lattice, we have first estimated the percolation 
threshold using the method described in Q . Estimates of 
p c in the thermodynamic limit are quite well-established 
for DP; however, we have worked with our estimated 
p c whenever finite size effects are strong, especially in 
three dimensions and rectangular geometries in two di- 
mensions. The shape of the lattice (i.e, when the size 
of the lattice is not equal in different dimensions) itself 
may also play an important role in determining, e.g., the 
spanning probabilities etc. Even in OP, it was shown || 
that the TB-spanning probabilities show scaling behav- 
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ior with respect to the aspect ratio. However, the perco- 
lation threshold is not affected except showing stronger 
finite size effects. With this in mind, and angular depen- 
dence being special to DP several variations have 
been introduced in two dimensions, e.g, two different di- 
rections of propagation have been considered and also 
the aspect ratio has been varied. In three dimensions, 
we only consider the DP to be propagating along the 
diagonal. 

Our main results can be summarized as follows: 

(1) We measure P(n) the probability for n TB-spanning 
clusters. From the absence of any significant finite size 
effect, we conclude that there is a nonzero probability of 
a finite number of such clusters at p — p c both in two 
and three dimensions. 

(2) The probability distribution P{n) ~ exp(— an 2 ) and 
the functional form is perhaps universal. 

(3) We conjecture that a is a weakly universal quan- 
tity that vanishes above the upper critical dimension and 
a = 0(e) for e = d c — d — > 0. This a is a new character- 
istic for percolation. 

(4) The fractal dimension of the spanning clusters is the 
same as expected for the unique cluster at p c . 

(5) The probability distribution of the masses (from all 
configurations) has a universal asymmetric Pearson dis- 
tribution (Type III) in terms of a scaled mass variable. 
See Eq. |. 

(6) The mass ratio of the two largest clusters also has a 
universal probability distribution. 

(7) The ratio of the average masses, universal for OP, 
may have an angular dependence. 

In two dimensions, we consider lattices of dimensions 
L x x L y . With L x = L y (square lattices), we check that 
for both directions, there is a non- vanishing probabil- 
ity for the existence of more than one spanning cluster. 
In the case of the diagonal DP, the number of spanning 
clusters n < 2 always, and the normalized probability 
P(2) = jj-pTo) here is found to be not significantly dif- 
ferent from that of OP, where the occurrence of coexisting 
spanning clusters is quite well established. This is shown 
in Fig |j]a and to be noted is the absence of any significant 
finite size effect. 

For the DP parallel to the y axis, a larger number 
of spanning clusters is obtained, giving us an oppor- 
tunity to study the behavior of P(n). Since in OP, 
P(n) ~ exp(— an 1 ) we check whether this behavior is 
still true for DP. We plot, in Fig ^, P{n) against n and 
n 2 and the latter gives a better straight line, at least for 
the range of n in the present study. Fig ^c shows a plot of 
the finite log-derivative defined as 7„ = ln(z) / ln(n/n+l) 
where z — ln(P(n+l))/ ln(P(n)) for two cases with n > 5 
for both d = 2 and 3. This plot of 7„ shown against 
though not very systematic, is indicative of an extrap- 
olated value of 2 for n — > oo and definitely rules out a 
small value like 7 = 1. That we obtain a larger number 



of clusters, compared to OP, may be due to a lesser value 
of a in DP. Its significance is discussed later. 

For rectangular geometries, i.e., when L x 7^ L yi we 
find that for the diagonal DP, P{n > 1) vanishes for 
p = L y /L x > 1. The plot P(n > 1) against p in Fig 
[j] for different sets of L x , L y shows that it is only de- 
pendent on the aspect ratio p, a result which was also 
obtained in OP l|. A larger number of spanning clus- 
ters is obtained even for the diagonal DP when p < 1. 
No such crossover is observed for the DP parallel to the 
spanning direction in the sense P(n > 1) is always non- 
zero even for p as high as 8. These results tell us that 
such a crossover behavior is dependent on the direction 
of the DP, apparently vanishing when it is parallel to the 
spanning direction. In fact, the crossover in the diagonal 
(and not in the parallel) DP is possibly because of the 
relatively larger spread of the clusters in the x direction. 
The study of the probability distribution P(n) is possi- 
ble for the diagonal DP with p < 1 as we get around 
5-6 spanning clusters, as well as for the parallel DP for 
any aspect ratio. We again get results compatible to the 
value 7 = 2 for each individual case. 

Another interesting behavior is observed for the DP 
propagating along the diagonal direction when we con- 
sider p < 1. Here p c decreases as L x is made larger, 
keeping L y constant. With L x /L y — m an integer, one 
may visualize the lattice as comprising of m square cells 
placed adjacent to each other. We see that though the 
value of p c is higher in each individual cell, it is lesser for 
the entire system. This could happen if the clusters in 
the individual cells do not TB-span but join with those 
in the neighboring cells to span the larger lattice. If this 
picture is correct then the width of the spanning cluster 
along the x direction should increase with L x . We have 
verified that this is true. 

The mass of the spanning cluster in a DP is expected 
to behave as (Mi^ n ) ~ L» 11 where Mi^ n is the mass of 
the ith largest spanning cluster with n spanning clusters 
present, L\\ is the length parallel to the preferred direc- 
tion, D\\ = —f3/v\\ + 1 + (d — l)/b, b = v\\/v±, v\\ (v±) 
is the correlation length exponents parallel (perpendicu- 
lar) to the direction of the DP, and (3 is the DP order 
parameter exponent Jj],0. -D|| — 1.47 in two dimensions 
with the values of (3, v\\ and v±_ as given in ]lq] . In the 
L x L lattice with the direction of propagation along the 
diagonal, it is expected that M^ n ~ i II as the spanning 
is considered from top to bottom. This is because the 
length scale in the direction of propagation is ~ L\\ and 
the width of the cluster is ~ L± and here Lu ~ L. We 
indeed find that D ~ 1.5 (comparable to 1.47) here. See 
Fig. |. For the DP parallel to the spanning direction, 
iii ~ L obviously, and again the above scaling behavior 
is observed. Only exception is the diagonal DP case with 
p > 1 where the mass exponent is found to take a new 
value ~ 1.8. A value close to 2 suggests a nearly compact 
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structure which could happen from strong boundary ef- 
fect in this particular geometry. The data of Fig || show 
that (Mi <n ) is independent of the number n of clusters 
present, signifying that the largest cluster is uncorrelated 
to the remainder of the lattice. 

The ratios of the masses of the individual spanning 
clusters should be constants as the scaling behaviors are 
identical. Defining r\ n ^ — (M^ n ) / ' (Mj >n ) , we find that 



is 1.36 ± 0.03 for the diagonal DP and 1.31 ± 0.02 



for the parallel DP for p = 1. The variations of these 
ratios are lesser when the aspect ratio is changed but the 
direction of the DP remains same. For comaprison, we 
add that rfl ~ 1.4 for OP @. 

In three dimensional symmetric lattices (£ 3 ), we find 
a significantly different larger number of spanning clus- 
ters compared to OP with the same number of ran- 
dom initial configurations (10 6 ). However, when we plot 
P(n) against n and n 2 , again the latter gives a better 
straight line. The extrapolated value of -f n is also pos- 
sibly close to 2 (see Fig ||). The fractal dimension of 
the cluster from the known values is around 1.67 and 
the log-log plot of Mij against L is compatible with 
this value (Fig || We obtain several mass ratios, for 
example, r{ 2 2 } = 1.52 ± 0.02, rf 2 } = 1.34 ± 0.03, and 

{3} 



r\ 3 = 1.77 ±0.05. As n becomes larger, the fluctuations 
also increase and the values of the ratios for n > 4 may 
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be less reliable. It maybe mentioned here that r\ 
for OP. Though our results are not inconsistent with the 
expectation of angular dependence in DP [|l(|, still it is 
perhaps premature to comment quantitatively on the uni- 
versality of these ratios from the present data. However, 
one can safely conclude that at least for three dimensions, 
these ratios are not same for OP and DP. 

In order to probe the detailed mass distribution, we 
obtain histogram for the cluster masses when DP is par- 
allel to y axis (the TB-spanning direction). A collapse 
of all the probability distributions is found, for vari- 
ous geometries and aspect ratios, if a scaled variable 
X = (Mi - (Mi)) /Si is used where Sf = (M?) - (M,) 2 , 
and Mi is the ith largest cluster irrespective of the value 
of n. No sensitivity to n was detected. We see a definite 
lower cutoff for the scaled variable. The plot is shown 
in Fig f| which also shows a Pearson distribution (type 
III). Since by construction the scaled variable X has a 
zero mean and unit standard deviation, only parameter 
available for the Pearson distribution is the lower limit. 
The resulting curve for — a < x < oo, 



Q(x) = 



r(« 2 ) 



[a(x + a)] c 



1 — a(x+a) 



(1) 



with a = 2.8, fits remarkably the probability distribution 
obtained numerically. 

The single variable distribution Q(x) does not reflect 
the restriction Mi > M^- For a joint distribution, we 
study the probability density lZ(y) for the mass ratio 



y = M\jMi of the two largest clusters, again irrespec- 
tive of n > 2. A data collapse with some finite size 
effect is seen, see Fig |]b. An exponential distribution 
TZ(y) = Aexp[— A(x— 1)] with A = 4. is also shown there 
for comparison. With this exponential distribution, the 
average mass ratio is expected to be 1 + A^ 1 » 1.25. 

Our main results have already been summarized at the 
beginning. We like to add a few comments here. First, 
the occurrence of the two simple distributions in Fig. |] 
remains unexplained. The second is regarding a in the 
distribution P(n) ~ exp(— an 2 ). As mentioned earlier, 
the number of TB-spanning clusters at the threshold is 
infinity at the upper critical dimension (d c = 5 for DP 
and 6 for OP), validity for all d < 5 for DP also. From 
P{n), Generalizing our result to all d (as known to be 
the case for OP), we may take n c = 1/a/ck as a measure 
of the number of clusters. This n c , defined strictly at 
the critical point, has to diverge as d — > d c . Therefore, 
a must decrease || with dimension, vanishing at d = d c . 
We indeed find that a for DP is lesser than OP (for three 
dimensions, where we could make a comparison), a fact 
compatible with d c (DP)< tf c (OP). Assuming analyticity 
in e = d c — d, one expects a = at + be 2 + .., for small 
e, with geometry (and OP/DP) dependent a and b. In 
fact, for OP, a = 0.15 and b — 0.05 give very good agree- 
ment with the a values obtained for e = 1,2 and 3 ||. 
We conjecture that critical fluctuations are responsible 
for nonvanishing of a, and therefore a is a new weakly 
universal critical quantity. 
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FIG. 2. (a) The probability distribution of spanning clus- 
ters plotted against n 2 for different DP's. (b) Same as (a) but 
now plotted against n. (c) The plot of j n vs 1/n. 
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FIG. 3. The masses of different spanning clusters for DP 
are shown along with specific power law variations for com- 
parison. For two dimensions, L — L x and system sizes are 
enlarged by 10 for three dimensions. 
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FIG. 1. (a) The normalized probability of the existence 
of two spanning clusters for ordinary and directed percolation 
are shown for square lattices. The direction of the DP is along 
the diagonal and P(n > 2) = for both cases. (b)The nor- 
malized spanning probabilities for n > 2 for the the diagonal 
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FIG. 4. (a) The probability distribution of masses of indi- 
vidual clusters irrespective of the total number. The solid line 
is the Pearson distribution type III with zero mean and unit 
standard deviation with a lower cutoff —2.8. (b)The probabil- 
ity distribution for the mass ratio of the two largest clusters 
for n > 2. The solid line is an exponential distribution with 
A = A. 
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